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We study the dissipative properties of a harmonic oscillator subject to two independent heat baths, 
one of which couples to its position and the other one to its momentum. This model describes a large 
spin impurity in a ferromagnet. We find that some effects of the two heat baths partially cancel each 
other. Most notably, oscillations may remain underdamped for arbitrarily strong coupling. This 
effect is a direct consequence of the mutually conjugate character of position and momentum. For 
a single dissipative bath coupled linearly to both position and momentum, no underdamped regime 
is possible for strong coupling. The dynamics of purity loss for one and two wave packets is also 
investigated. 
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The harmonic oscillator provides a scenario where the 
physics of a quantum particle coupled to a heat bath 
can be investigated analytically 0, 0, • Thus its study 
may shed light on physical effects that have been found 
in less tractable models. Recently, Castro Neto et al. 
4] have investigated, by means of perturbative renoma- 
lization group methods, the equilibrium dynamics of a 
quantum magnetic impurity in a ferromagnet and have 
found that the symmetric coupling of the magnon bath to 
the x and y spin components (z being the magnetization 
direction) effectively decreases the strength of the cou- 
pling. They have coined this term quantum frustration, 
since it may be interpreted as the inability of the spin 
bath to simultaneously measure two non-commuting ob- 
servables such as the x and y components of the impurity 
spin. An important consequence is that spin coherence 
may be longer lived than it would if only one of the spin 
components were coupled to the dissipative bath. This 
effect could be relevant for quantum information, as it 
would provide a mechanism for quantum spins to remain 
coherent for long times. 

Here we investigate the effect of quantum frustration in 
a different but exactly tractable physical system, namely, 
a harmonically oscillator coupled separately, through its 
position and momentum, to two independent oscillator 
heat baths In the symmetric limit, this model de- 
scribes the behavior of a large spin magnetic impurity 
in a ferromagnet. Another instance of a physical sys- 
tem which interacts, through position and momentum, 
with two different baths is a Josephson junction, where 
the relative particle number (proportional to the electric 
dipole) couples to the radiation field while the relative 
phase interacts with the quasiparticle field 0. Those 
two environments have different spectral properties. To 
investigate a more symmetric coupling, we consider a har- 



monic oscillator whose position and momentum interact 
linearly with two independent Ohmic baths Q. We find 
that the two baths do cancel in some but no all respects. 
A degree of cancellation is revealed by the persistence of 
underdamped oscillations for arbitrarily strong dissipa- 
tion provided that the two baths couple with comparable 
strength. Quantum purity is weakened by the presence 
of a symmetric second bath although its decay is slowed 
down. Since a symmetrically damped harmonic oscilla- 
tor behaves like a large spin in a ferromagnet, we refer 
to such a mixed situation as quasiclassical frustration. 
We investigate the following model Hamiltonian Q: 



H = V q (q+Sq)+V p (p+Sp} 
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where [q,p] = i and g q ,g p are sufficiently well behaved 
functions. For general g q ,g p the two equations in Eq. PJl 
cannot be decoupled without generating interactions be- 
tween the two baths. However, if one of the coupling 
functions, say g p , equals 1, then [Sp, Sq] — and it be- 
comes possible to remove both fluctuating contributions 
from the potential terms through the unitary transfor- 
mations Up = enp(ipSq) and U q — exp(i J q+Sq Sp(q')dq'). 
One arrives at the Hamiltonian 

H = V q (q) + J2^qk 



+ V P (p) 
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where the short-hand notation \a\ 2 = a 1 a has been used. 

At this point it is important to specify what we mean 
by coupling to position or momentum. Those are the 
particle variables to which the environment couples as a 
set of otherwise independent harmonic oscillators. This 
is the case e.g. in ©, where the bath oscillators would 
remain independent if q and p were c-numbers. On the 
contrary, this is not the case in (JJJ, where, due to the 
nonlinear character of 14_and V p , the bath oscillators do 
interact with each other ja] ■ A popular example is that of 
a charged particle interacting with the photon field [6j. 
When the velocity- coupling model is adopted, it must 
be accompanied by a diamagnetic term that contains 
an interaction between photons. By contrast, through a 
canonical transformation, the photon field interacts with 
the particle position without an explicit intcrphoton in- 
teraction. Thus, within the precise convention we pro- 
pose here, the electromagnetic field couples to the posi- 
tion of a charged particle. 

We focus on the case of a harmonic oscillator [V q (q) = 
u> q q 2 /2, V p {p) = co p p 2 /2} with g q = g p = 1: 
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This is the model of dissipative coupling that displays 
the highest degree of symmetry between q and p. Hence 
it is interesting to study frustration in a physical system 
other than a quantum spin in a ferromagnet. The two 
baths of independent harmonic oscillators are described 
by the spectral densities 



J n {w) = 2^\C nk \ 2 S(uj-uj nk ) 



p,q 



(5) 



We assume a power law behavior at to — and write 
J„(cj)= 27„(jj Q "/(tjp^~ 1 7r), where the introduction of the 
frequency ui p h renders the coupling constants 7„ dimcn- 
sionless. Moreover, large cutoff frequencies Q„ are as- 
sumed to exist for both environments. 

Eliminating the bath variables, the Heisenberg equa- 
tions of motion for q and p are obtained, 



q(t) = U p p(t) + 



J dsK p (t 



s)p(s)+F p (t) 



-**)=«,?(*)+ J dsK q {t- S )q(s)+F q (t), (6) 

where K n {t) = J n (uj) cos(ujt)d\iiuj and F n (t) = 
Y^,k CnkCtnk exp(-iw n fet) + H.c. In Fourier space, Eq. © 
reads 



J q (tL>) - LOq 



itop = F q 



—iwq + J p (u>) — ojp p = F p 
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where J n (uj) is the symmetrized Riemann transform jtj: 



y ^-zsgnH^/(|c|) 



(9) 

The oscillation modes are given by the zeros of the 
polynomial 

X _1 M = Ujl-UJ 2 -UJ q J p (u!)-U!pJ q {u!) + J q (uj)Jp(uj) (10) 

where x( w ) is the generalized susceptibility. 

We further assume that the two heat baths are^Ohmic: 
J n (io) = 2"/ n oj /tt. For Q n — ► oo, this implies J n (u>) — > 
ij n io. Then the eigenfrequencies are given by 

0Jq - i{u q lp + w p 7 9 )w - (1 + 7g7 P )w 2 = , (11) 

the solutions being 
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The transition from k < 1 to n > 1 marks the crossover 
from undcrdamped to overdamped oscillations. The con- 
dition k < 1 requires (criterion A) 



\lqUp ~ J q UJ p \ < 2u 



(14) 



The underdamped region satisfying (|14|) lies in a stripe of 

width A = Ar\ (l + ?/ 4 ) 1 ^ 2 with r\ — (ojq/ujp) 1 ^ 2 , limited 
by the graphs of the functions f(j q ) = (7 g ±2?y)f7~ 2 . The 
stripes of underdamped oscillations in the (7 9 ,7 P ) plane 
are plotted in Fig. ^ for 77 = 1/3, 1, 3. A remarkable con- 
sequence is that, given a value of e.g. j q , one may drive 
the system from the overdamped to the underdamped 
regime by increasing 7 P . For 77 = 1, the oscillator is un- 
derdamped if j q ~ 7 P , i.e., if the couplings to the two 
baths are of comparable strength. When "f q = 7 P = 7 
the dimensionless parameter k becomes, 



7 (l + 7 2 ) 
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Thus, in the fully symmetric case, the oscillator remains 
underdamped for all values of the coupling strength. 
This is in contrast with the case of one noise [pf p = 0), 
which requires j q < 2t) to be underdamped. Conversely, 
if -f q = 0, the condition is j p < 2r\~ x . The inset of Fig. 
shows the underdamped region for an oscillator coupling 
through q and p to a single heat bath ^(J . Surprisingly, 
the behavior is qualitatively different in that underdamp- 
ing is always lost for sufficiently large coupling. 

Another interesting quantity is D q (uj) = Imx gg (u;)/a;, 
where Xqqfa) is the Fourier transform of ([q(t) , q(0)]) . 
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FIG. 1: Stripes of underdamped oscillations in the (7 9 ,7 P ) 
plane for three different values of the parameter 77 = 
(uq/tjp) 1 '' 2 , for q and p coupled to different baths. Inset: 
Shaded region denotes underdamping for a single bath cou- 
pled to q and p. 



For Ohmic environments, 



D q (u) 



Iqul +7 P (1 + 7g7 P )^ 2 



[(1 + Iqlp)^ 2 ~ W 2 ] 2 + (iqUp + JpUq) 2 U) 2 
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Following Ref. we may view the presence of a peak 
in D q (uj) as a signature for the existence of coherent dy- 
namics (criterion B). This occurs for 



(17) 



In the symmetric case, this translates into 7 < \/3. For 
an oscillator coupled to a single bath through its position 
(lp = 0)> the requirement is j q < \/2rj. Finally, if 7 9 = 0, 
D q (ui) always displays a peak |ll|. 

A third possible condition for the existence of coherent 
dynamics is that, in Eq. 1|12|) . 



I Im w± I < I Re to± | . 

(criterion C). This yields 

22, 2 2 „ r, 2 

IqUp+JpUq < 2LU 



(18) 



(19) 



In the symmetric case, the condition l|19|) becomes 7 < 1. 
For j p = 0, it becomes j q < \f%], while for 7g = 0, it 
reads 7p < V^i] 1 ■ 

In table [I] the coherence signatures for the main three 
particular cases are summarized. Criteria A and C are 
symmetric in q and p, but not criterion B. The general 
trend (particularly clear if one considers A and B) is that, 
starting from a single dissipative bath coupled to e.g. 
q, the introduction of a second bath that couples to p 
with the same spectrum and comparable strength favors 
coherent and underdamped dynamics. For example, an 





symmetric 


only 7, 


only 7 P 


K < 1 (A) 


always 


7s < 2?7 


7p < 27T 1 


D q (cj) has peak (B) 


7 < V3 


7g < V%V 


always 


\Tmu±\ < |Rew±| (C) 


7 < 1 




7 P < a/27/- 1 



TABLE I: Condition for coherent dynamics according to three 
different criteria (left column, labeled A, B, C), for three par- 
ticular cases (upper row). The symmetric limit includes the 
assumption 77 = 1. General case is given in main text. 



oscillator with 77 = 1 that is driven from j p — to 7 P = 
7q = 7, with j q fixed at a value \f2 < 7g < \/3, will 
cross over from incoherent to coherent behavior under 
both criteria A and B. 

We have noted the striking result that, in the sym- 
metric case, the oscillator is underdamped for all values 
of 7. However, criteria B and C indicate that the os- 
cillator lies deep in the underdamped region only if 7 
is small. Such a limitation is also patent in the (pos- 
sible) maximum of D q (u>) / D q (Q) as well as in the ratio 
I Rew±|/| Imw±| = 7 _1 . Both quantities stay well above 
unity only if 7 is small. A related point is that, as 7 — > 00, 
the ratio D q (oj)/ D q (0) does not saturate but rather de- 
cays as Wg/7 2 o; 2 for nonzero u>. 

An oscillator initially prepared in a pure coherent state 
that at t = begins to interact linearly with an oscillator 
bath is described by a reduced density p which remains 
Gaussian at all times. Then the purity V(t) = Tr(p 2 ) is 
given by 

V- 2 (t)=4{q 2 ){p 2 ) . (20) 

At long times, (q 2 ) and (p 2 ) reach their equilibrium 
values, which contain contributions from both baths as 
well as hybrid terms which vanish if any of the two baths 
disappears For simplicity, we focus on the zero tem- 
perature, weak coupling case. Then, 

I 7<J A fin 1" 



•7p 



ln^ 

UJ 



■0(7„7p), ( 21 ) 
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and similarly for (p 2 ). The two baths have opposite ef- 
fects on (q 2 ), but for Q n 3> ujq the logarithmic divergence 
from the p-coupled bath overwhelms the squeezing of q 
favored by the g-coupled environment. An impure mix- 
ture results: 

2j q ( '. f2„ 



1 + zj* I in 12 
2j p ?i (in 
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LO 



1 +0( 7g>7p )»l. (22) 



The dynamical evolution of the purity V(t) is cumber- 
some in the general case, but it becomes tractable in the 
symmetric problem. We find (f2 g = fl p = tt) 

( „-fii 
27 



T(t) 



o<t<n- 



1 + 



(l +7 2)3/2 



Q~ L < t 



00 



(23) 
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with = 70)0/(1 + 7 2 ) = |Imo>±|. The fast decay on 
the scale of fi -1 comes from the choice of decoupled ini- 
tial conditions 01 an d horn the concurrence of two baths 
[l3|. For t the system evolves slowly towards 

equilibrium. The divergence of t for 7 — > 00 might be 
interpreted as robustness against purity loss. However, it 
should be noted that such a slowing down merely reflects 
a dilation of all time scales with increased friction. For 
instance, |Reo±| = o /(l + 7 2 ) vanishes even faster. 

We have also investigated the purity decay when at 
t = the system is prepared in a linear superposition 
of two coherent states centered at q — ±a/2 with zero 
average momentum 14]. For the symmetric case, we find 

M 



p(t) 



Pit) 



cosh 



cosh^ (a 2 /8) 



(24) 



where <p(t) is a complicated function that evolves from 
4>(0) = 1 to limt^oo <p(t) = 0, and the single wave packet 
purity V(t) is given in As expected, P(t)/V(t) -> 1 
as a -> 0, and P(t)/V{t) -> 1/2 as a -> 00. 

Interestingly, the structure of l|24|l is such that, as 
time passes and 4>{t)P(t) evolves from 1 to 0, the ra- 
tio P(t)/P{t) starts at unity, as corresponds to a pure 
state, then decreases and finally, at long times, goes back 
to unity. When a is large P(t)/P(t) decays rapidly on 
a timescale ~ l/4a 2 7 to 1/2. There it stays for a time 
which increases with distance as ~ 7 _1 In a. Afterwards 
it returns to one. The ratio 1/2 can be rightly inter- 
preted as resulting from the incoherent mixture of the 
two wave packets. Thus it comes as a relative surprise 
that P(t)/P(t) becomes unity again at long times, as 
if coherence among the two wave packets were eventu- 
ally recovered. The physical explanation lies in the er- 
godic character of the long time evolution, with both 
wave packets evolving towards the equilibrium configu- 
ration described in Eqs. 1201221 Once the two initially 
separate wave packets begin to overlap, they regain mu- 
tual coherence. Due to the symmetry of the problem, a 
similar result would have been obtained if the oscillator 
had started from a superposition of two coherent states 
located in the same region of real space but with different 
average values of the momentum. 

In summary, we have found features that are remi- 
niscent of an effective particle-bath decoupling, such as 
the persistence of underdamped oscillations for arbitrar- 
ily large values of 7 in the case of a symmetric oscilla- 
tor and the slowing down of purity decay for a Gaussian 
wavepacket. Another feature is that two initially sepa- 
rate wave packets regain relative coherence at long times 
because they recombine. The situation is reminiscent of 
the quantum frustration exhibited by a magnetic quan- 
tum impurity albeit in a more limited form 41. The mam 



difference between the two problems is the dimensionality 
of the particle Hilbert space. Evolving in the continuum, 
the quantum oscillator can be considerably degraded by 
the effect of the environment, as shown in 1)23(1 . It is only 
when the two wave packets recombine because of ergod- 
icity that relative purity is recovered. By contrast, the 
spin-1/2 magnetic impurity lives in a two-dimensional 
space. The only possible effect of the environment is to 
flip the spin. Thus, in any representation, two initially 
orthogonal states quickly overlap and tend to preserve 
mutual coherence. The net result is an increased decou- 
pling from a symmetrically dissipative environment. The 
requirement of low dimensionality suggests that strong 
frustration is a genuinely quantum effect. 
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